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Section  1 
INTRODUCTION 

The  objective  of  this  study  is  to  determine 
if  a  basis  can  be  established  for  a  statistical  model  of 
the  distribution  of  vertical  shear  (vertical  gradient  of 
horizontal  velocity)  in  the  upper  ocean.  We  have  approached 
this  problem  in  three  stages: 

•  Define  as  our  basic  descriptor  the  probability 
distribution  of  shear  squared; 

•  Determine  how  well  this  observed  probability 
distribution  can  be  described  by  a  known 
analytic  probability  distribution; 

•  Test  the  assumptions  which  would  have  to 
be  satisfied  to  ensure  that  the  observed 
samples  do,,  indeed,  come  from  the  known 
distribution. 

In  order  to  examine  these  questions,  we  have 
analyzed  vertical  profiles  of  horizontal  velocity  in  the 
upper  ocean  provided  by  David  Evans  of  the  University  of 
Rhode  Island.  The  data  were  collected  by  a  sensor  system 
called  YVETTE. 

In  Section  2  we  briefly  describe  the  YVETTE  data. 
A  more  detailed  description  is  given  by  Lambert  et  al . 
(1980).  We  also  describe  the  processing  which  we  applied 
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to  the  YVETTE  profiles  and  the  general  characteristics  of 
the  data.  In  Section  3  we  derive  an  analytic  model  for  the 
probability  distribution  of  mean-square  shear.  This  model 
leads  us  to  speculate  that  the  occurrence  of  shear-squared 
values  in  various  depth  regimes  follows  ax2  probability 
distribution.  In  Section  4  we  test  this  hypothesis 
statistically  in  two  ways.  In  one  we  compare  directly  the 
observed  shear-squared  distributions  with  suitable  x  2- 
distnbutions  to  assess  the  acceptability  of  the  hypothesis, 
while  in  the  other  we  test  how  well  the  statistical  (as 
opposed  to  physical)  assumptions  required  for  a  x ^distri¬ 
bution  are  satisfied  by  the  data.  In  Section  5  we  describe 
an  attempt  to  implement  the  x2  model  by  making  use  of 
a  possible  relationship  between  n2  and  S2  proposed  by 
Patterson  et  al.(1981)  in  a  companion  SAI  report.  We 
discuss  the  utility  of  this  model  and  make  recommendations 
for  further  improvements. 


Section  2 


YVETTE  VELOCITY  PROFILES  -  OVERVIEW 

2.1  YVETTE  STATIONS 

The  data  on  which  our  analysis  is  based  were 
obtained  by  David  Evans  of  the  University  of  Rhode  Island  in 
1975  and  1977  with  a  vertical  profiling  instrument,  YVETTE. 
At  the  time  these  profiles  were  measured  YVETTE  was  a 
four -meter-1 ong  tube  equipped  with  Neil  Brown  conductiv¬ 
ity,  temperature,  and  depth  (CTD)  sensors  as  well  as  an 
orthogonal  pair  of  acoustic  current  meters.  A  detailed 
description  of  YVETTE  is  provided  by  Evans  et  al .  (1979). 

When  deployed,  the  probe  fell  freely  through  the  water 
column  at  a  rate  of  25  cm  s~l.  The  water  temperature, 
conductivity,  two  orthogonal  components  of  horizontal 
velocity  relative  to  the  probe  housing,  and  the  orientation 
of  the  instrument  relative  to  the  earth's  magnetic  field 
were  recorded  internally  on  magnetic  cassette  tape.  The 
sensors  were  sampled  at  a  rate  of  2.5  times  per  second  to 
yield  a  vertical  sampling  interval  of  approximately 
10  cm.  The  CTD  sensors  had  a  resolution  of  0.001  mmho 
cm-1,  0.0005°C,  and  0.05  dbar,  respectively.  The 

acoustic  current  meters  had  a  sensitivity  of  0.05  cm  s_1 
and  the  compass  could  resolve  orientation  to  within  3°. 

The  data  set  provided  to  us  was  filtered  to  remove 
variations  on  vertical  scales  less  than  2  m.  The  velocity 
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profiles  were  corrected  L;  horizontal  drag  effects  imposed 
on  the  system  housing  by  the  large-scale  vertical  shear. 
A  description  of  this  processing  is  contained  in  Lambert 
et  al.  (1980).  The  resulting  data  set  should,  in  principle, 
include  velocity  structure  resolved  over  vertical  scales 
between  2  m  and  100  m. 

The  data  set  consists  of  fifteen  stations, 
fourteen  in  various  parts  of  the  North  Atlantic  Ocean  and 
one  in  the  Pacific  Ocean.  The  position,  date,  and  time  of 
these  profiles  are  listed  in  Table  2.1.  Figure  2.1  shows 
the  profile  locations. 

2.2  VELOCITY  DATA  PROCESSING 

The  profiles  of  horizontal  velocity  obtained  from 
Evans  consisted  of  values  spaced  irregularly  in  the  vertical 
(pressure)  direction  with  a  typical  spacing  of  one  meter. 
Our  first  step  in  analyzing  shear  values  was  to  inter¬ 
polate  the  irregularly  spaced  values  to  a  profile  of  values 
at  regular  intervals.  We  did  this  with  a  cubic-spline 
interpolation  to  produce  values  of  u  and  v  at  intervals  of 
lm.  We  then  performed  numerical  differentiations  over 
specified  separations.  We  considered  separations,  Az,  of  2, 
4,  8  and  16  m  in  this  study. 

We  define  shear-squared  at  depth  z  as  the  squared 
magnitude  of  the  difference  between  the  horizontal  velocity 
vector  v  at  depth  z  +  Jaz  and  z  -JAZ,  normalized  by  (Az)2. 
That  is 


S2(z;az)  =  |  v(z+^Az)  -  v(z-iAz)  |2(Az)“2  , 

-  (Azu)2  +  (Azv)2  ,  (2.1) 
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»  Table  2.1 

YVETTE  STATIONS1 


Station 

Number 

Time 

(GMT) 

Date 

Latitude 

(N) 

Longitude 

(W) 

Ccranent 

5 

1942 

5  Nov.  75 

32°19' 

64° 34' 

Near  Bermuda 

8 

0225 

8  Nov.  75 

35°00' 

66° 30' 

Sargasso  Sea 

9 

1219 

It 

II 

It 

II 

10 

1814 

9  Nov.  75 

38°09' 

69°06' 

Gulf  Stream 

11 

0036 

10  Nov.  75 

38°05’ 

69°03 ' 

If 

12 

1312 

II 

38°15' 

69°07 ' 

If 

18 

— 

7  May  77 

22°47 1 

70°43 ' 

Edge  of  thermocline  eddy 

21 

— 

9  May  77 

22° 27' 

70° 57' 

Center  of  thermocline  eddy 

23 

— 

16  May  77 

36°24' 

67°36' 

Outer  part  of  GSR^ 

24 

— 

II 

36°20' 

67°  44' 

Midway  along  radius  of  GSR 

25 

— 

17  May  77 

36°09' 

67°53' 

Near  center  of  GSR 

NOR  1 

— 

1973 

— 

— 

Norwegian  Fjord 

NCR  4 

— 

II 

— 

— 

II 

NCR  6 

— 

11 

— 

— 

Norwegian  Coastal  Current 

EP0CS36 

— 

July  79 

0°03' 

109°57 ' 

Data  start  just  below  EUC4 

1  Adapted  from  Lambert  et  al.  (1980) 

2  Gulf  Stream  Ring 

3  Equatorial  Pacific  Ocean  Climate  Study 

4  Equatorial  Undercurrent 
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Figure  2.1  Positions  of  YVETTE  Stations 


where  u  and  v  refer  to  velocity  components  in  the  east-west 
and  north-south  directions  and 


l2u  =  |u(z^iz'<- 

with  a  similar  expression 
overlapping  differences,  we 
as  we  increased  the  spacing. 


u (.z-jAz  iz  ^  (2.2) 

for  Azv .  Since  we  used  non- 
obtained  fewer  shear  estimates 


Once  the  Azu  and  A  zv  values  were  calculated,  we 
calculated  the  depth  average  of  each  profile,  defined  as 


/ 

“  z 


u 


1 

K 


I 

i=l 


1zu(zi) 


(2.3) 


with  a  similar  expression  for  Azv .  By  subtracting  the 

average  from  the  individual  values  we  obtained  modified  data 

A  u ’  and  A  v * : 
z  z 


Azu'(zi)  =  Azu(zi)  -  AZU  .  (2.4) 

In  this  way  we  removed  any  linear  trend  which  might  have 
occurred  in  the  data. 

The  S2  values  used  in  this  study  are  defined  as 
S2  =  (Azu’)2  +  (Azv’)2  .  (2.5) 
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SAMPLE  STATION  -  YVETTE  08 


The  first  station  we  analyzed  was  YVETTE  08 
from  the  Sargasso  Sea.  Figures  2.2  and  2.3  show  plots 
of  the  data  received  from  URI  for  this  station.*  Note  the 
presence  of  an  upper  isopycnal  layer,  a  sharp  pycnocline 
below  that,  and  a  region  of  less  steep  but  still  increasing 
density  below  that.  We  have  divided  the  profiles  for  this 
station  into  three  regimes:  the  surface  mixed  layer  (ML), 
the  upper  tbermocline  (UT),  and  the  deep  layer  (DL).  The 
dotted  lines  in  the  figures  denote  rather  subjective  bound¬ 
aries  between  these  layers.  The  values  of  A  2u  and  A2v 
tend  to  be  highest  in  the  upper-thermocl ine  region  where  the 
Brunt-Vaisala  frequency  is  highest  and  lower  in  the  surface 
mixed  and  deep  layers.  This  apparent  correlation  is  ex¬ 
plored  in  detail  by  Patterson  et  al.  (1981).  We  will  return 
to  its  implication  later  in  this  document. 

All  of  the  YVETTE  profiles  were  divided  into  sub¬ 
profiles  in  this  manner,  using  divisions  established  by 
Patterson  et  al .  (1981)  on  the  basis  of  changes  in  Brunt- 
V&is'Aia  frequency.  We  have  retained  their  definitions 
of  the  layers  without  modification. 


*  The  profiles  shown  in  these  figures  were  taken  directly 
from  data  provided  by  Evans.  The  actual  data  points  are 
located  at  irregularly  spaced  depths  as  mentioned  earlier. 


(a)  East  Velocity  Component 


(b)  North  Velocity  Component 


(c)  Current  Speed 


(d)  Current  Direction 


(e)  Azu 


( f  )  A2v 


Figure  2.3 


Profiles  of  the  (a)  east  and  (b)  north  component 
of  velocity,  the  (c)  magnitude  and  (d)  direction 
of  the  horizontal  velocity  vector,  and  the 
vertical  gradient  of  the  (e)  east  and  (f)  north 
component  of  velocity  for  YVETTE  station  8  in 
the  Sargasso  Sea.  The  dotted  lines  are  at  135 
and  250  dears. 


Shear-squared  profiles  which  we  computed  by 
applying  (2.5)  to  the  interpolated  velocity  profiles  from 
this  station  are  shown  in  Figure  2.4.  The  calculations  with 
2,  4,  8  and  16  m  spacings  are  all  presented. 

In  order  to  examine  the  distribution  of  values 
of  S2  we  divided  the  range  from  S2  =  o  to  S2  =  maximum 
in  each  profile  segment  into  equi-width  bins  such  that  the 
range  would  be  subdivided  into  30-40  bins.  We  then  gene¬ 
rated  histograms  of  the  number  of  values  falling  into  each 
bin.  These  were  expressed  as  the  fraction  of  the  total 
number  of  v?  acs  in  each  profile  segment.  Histograms  for 
the  three  depth  regimes  of  YVETTE  08  are  shown  in  Figure 
2.5a-c.  The  histograms  indicate  that  the  data  display  a 
strong  preference  for  low  values  of  S2  compared  to  the 
total  range  of  observed  values.  The  coarse  quality  of  the 
histograms  for  8  m  and  16  m  separations  is  due  to  the 
small  number  of  S2  values  in  the  samples. 
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Figure  2.4  Shear-squared  profiles  for  spacings, 
£z,  of  2,  4,  8  and  16  m  from  inter¬ 
polated  YVETTE  08  data. 
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SHEAR  DISTRIBUTION  MODEL 


It  is  possible  to  derive  a  model  for  a  shear 
probability  density  function  B(S2;az)  with  the  application 
of  a  few  assumptions  about  the  statistical  properties  of  the 
shear  field.  The  assumptions  are  simple  (but  not  trivial). 
In  this  section  we  derive  such  a  model  and  in  Section  4  we 
apply  it  to  the  data. 

The  basic  idea  was  developed  by  Bretherton  (1969). 
He  supposed  that  only  one  component  azu  of  the  vertical 
velocity  difference  was  significant  and  that  this  differ¬ 
ence  represented  the  sum  of  the  vertical  gradients  asso¬ 
ciated  with  a  large  number  of  linearly  independent,  random 
internal  waves.  Thus,  by  the  central -1 imit  theorem,  a  zu  is 
normally  distributed,  and' 

B(A2;U;Az)  =  (aV2  tt  )  -1  exp[  -  £(  Azu)2a-2]  (3.1a) 

where  the  variance 

o2  =  < (  Azu) 2>  (3.1b) 

and  the  braces  indicate  ensemble  averaging  (or  over  a  very 
large  number  of  independent  samples). 

The  density  function  is  defined  in  the  usual  way 

with 


Prob  (  Aq<X<1  0+d  X)  =  B(X0 ;  Az)  dx  . 


(3.2) 


The  probability  density  function  for  S2(Az)  (where  sub¬ 
scripts  1  refers  to  a  single  component  of  velocity)  is 
given  by  the  chi-square  density  with  one  degree  of  freedom: 

B(Si;Az)  =  [a  y/2~Sl  r(4  >  1  ^expC-iSia"2 ) .  *  (3.3) 

We  extend  the  analysis  to  allow  two  velocity 
components  to  contribute  to  the  mean-square  shear. 
We  assume  that  the  velocity  difference  vector  (Azu>  Azv) 
represents  the  vector  sum  of  velocity  differences  associated 
with  a  large  number  of  linearly  independent  random  phe¬ 
nomena.  Then  the  Multidimensional  Central  Limit  Theorem 
(Loeve,  1960)  implies  that  the  two-dimensional  proba¬ 
bility  density  function  B(azu,azv;Az)  of  the  vector 
approaches  a  two-dimensional  normal  distribution.  Thus, 


B(  A^u ,  A^v ;  Az)  =  (2ir 

.2 


z  '  z 
x  exp 


-i(l-r2)o  2 


(AZU) 


+  (Azv)‘ 


(3.4) 


where  a  ^  is  as  defined  in  (3.1b)  and  r  is  the  covariance 
r  =  <Azu  azv>  a-2 

We  assume  that  the  wave  field  is  isotropic. 

Therefore. 


r  =  0 


(3.5) 


*  r(  £)  is  approximately  1.772  so  that  the  quantity  in  the 
brackets  is  approximately  2.5a  S^. 
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and 


<AzuAzv>  =  0  • 


Each  individual  probability  density  is  given  by  (3.1)  so 
that  the  velocity  components  are  normally  distributed.  The 
separability  of  the  probability  densities  is  a  definition  of 
statistical  independence.  Hence  the  velocity  difference 
components  £zu  and  £zv  are  not  only  uncorrelated  (3.5) 
but  statistically  independent. 

It  is  now  straightforward  to  obtain  the  proba¬ 
bility  density  function  of  S2.  The  sum  of  the  squares  of 
n  normally  distributed,  independent  random  variables  with 
identical  distributions  has  a  chi-square  probability 
density,  with  n  degrees  of  freedom  (Papoulis,  1965).  Thus 
S2  has  a  2-degree-of-f reedom  density  given  by 


B(S2 ; Az)  =  (2o2 )  1  exp  (-iS2a  2) 


(3.6a) 


for  S2>  0.-  The  density  function  in  (3.6)  is  not  sharply 

peaked  for  small  values  of  its  argument  whereas  that 
in  (3.3)  is.  At  large  values  of  its  argument,  the  value 
given  by  (3.6a)  is  larger  than  given  by  (3.3). 

We  can  integrate  (3.6a)  to  show  that  the  expected 
value  of  S2  is 


^S2  =  <S2> 


(3.6b) 


=  2a‘ 


and  that  the  variance  of  S2  is 

o2  =  <S2  -  <S2>  >  = 
S2 


(3.6c) 


We  propose  that  (3.6a)  be  tested  as  a  model  for 
the  probability  density  of  mean  square  shear.  Instead  of 
B(S2;Az)  a  more  useful  statistical  description  is  given 
by  the  probability  distribution  function 


X2(S2;  Az)E 


;  Az)dS2 


(3.7) 


so  that 

Prob  (S2iSo2)  =  x2(S02;Az). 
From  (3.6)  and  (3.7) 


X2(S2;Az)  =  l-exp( -JS2o  2). 

(3.8) 

Figure  2.5  referred  to  earlier  shows  the  S2 
histogram  computed\  from  data  from  YVETTE  08  data.  Super¬ 
posed  on  that  figure  are  B(S2;  Az)  density  functions  given 
by  (3.6).  In  theory,  the  component  variance  used  in  the  x2 
probability  is  given  by 


a 


2  = 


2  2 

That  is,  ax  and  Oy  (assumed  identical)  are  the  variances 

of  the  population  distributions  of  the  shear  components. 

A 

In  practical  terms,  however,  we  used  the  estimate  a2  given 
by 


/N  /\ 

c2  =  i(a2 


'V 

°y  } 


(3.9) 


the  average  of  the  estimates  of  the  variances  of  the  two 
components  of  the  sample,  since  the  two  component  variances 
are  not  exactly  equal. 


t 


a 


4 


i 


/N  /\ 

2  2 

The  definitions  of  ax  and  a y  lead  to 


N  N 

—  57  (Au./Az)2  +  —  57  (Avi/A: 


(3.10) 


where  N  is  the  number  of  samples  in  the  depth  regime.  Note 
that  the  means  of  the  shear  components  have  been  removed 
(i.e. ,  iu.are  equivalent  to  Au^  of  2.4). 

Combining  (3.10)  and  (2.5)  leads  to: 

^  n  r  -i 

02  =  h  —  V*  [(Au./Az)2  +  (Av./Az)2J 


_  ,  N  /N 
*  N-l  yg2 


(3.11) 


Section  4 

SHEAR  DISTRIBUTION  ANALYSIS 

In  Section  3  we  derived  a  x^  model  for  the  dis¬ 
tribution  of  upper-ocean  mean-square  shear.  In  this  section 
we  test  the  applicability  of  that  model.  Our  approach  makes 
use  of  the  tools  of  hypothesis  testing,  in  particular,  a 
Kolmogorov-Smirnov  "goodness-of-f it"  test. 

We  characterize  the  appropriateness  of  the  model 
in  terms  of  the  critical  level  of  significance  (defined 
below)  ,  which  we  can  attach  to  the  hypothesis  that  the  S^ 
values  are  drawn  from  a  population  which  has  a  x2  distri¬ 
bution.  We  estimate  the  population  variance  as  the  variance 
of  the  sample.  A  goodness-of-f it  test  for  a  hypothetical 
distribution  in  which  some  significant  aspect  of  the 
hypothesized  distribution  is  calculated  from  the  sample 
is  a  "conservative"  test.  This  means  that  the  test  results 
are  biased  toward  higher  critical  levels  of  significance.* 

We  also  describe  the  results  of  (nonconservative) 
statistical  tests  performed  on  the  data  to  ascertain 
whether  the  conditions  on  the  velocity  field  necessary  for 
a  x2  shear  distribution  were  met.  We  applied  a  Cox-Stuart 
trend  test  for  randomness  of  the  samples  of  the  individual 


*  Conservative  here  implies  that  since  the  results  are 
biased  toward  higher  apparent  levels  of  significance  we 
are  less  likely  to  make  a  decision  to  reject  a  stated  null 
hypothesis. 
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components  Azu  and  AzV;  a  Lilliefors  test  for  normality 
of  the  components;  a  x2  test  for  mutual  independence  of 
the  components;  and  a  Smirnov  test  for  identical  distribu¬ 
tion  functions  of  the  two  components. 

A  review  of  hypothesis  testing  is  presented  in 
Appendix  A. 

4.1  GOODNESS-OF-FIT  TESTS  FOR  SHEAR 

4.1.1  Critical  Levels 

In  order  to  quantify  the  degree  to  which  shear 
distributions  for  each  profile  segment  can  be  represented 
by  a  x 2  distribution  (conservatively,  with  variance  set  to 
the  sample  variance),  we  applied  the  Kolmogorov  goodness-of- 
fit  test.  The  sole  assumption  required  by  this  test  is  that 
the  values  of  S2  be  random. 

The  test  consists  of  comparing  the  observed  prob¬ 
ability  distribution  G(S2),  calculated  from  the  sample,  to 
the  hypothesized  distribution  X2(S2)  with  a2  estimated 

A 

by  a2*  We  obtained  the  observed  probability  distribution 
simply  by  ordering  S2  values  from  low  to  high  (index  j)  and 
calculating  the  cumulative  sum  of  values  less  than  or  equal 
to  a  given  value: 


N 

G(S?)  =  i  H(xi  -  x  )  where  H«l,  Xj  <  xi 

j=l 

=0,  Xj  >  x.  .  (4.1) 

A 

o  °  2 

The  corresponding  hypothetical  probability  is  x2^i>a 
The  test  statistic  is  the  maximum  difference  between  these 
two  functions  at  any  Sit 


T 


sup 

«! 


G(sf) 


(4.2) 


where  "sup"  refers  to  the  maximum  value.  The  empirical  T 
is  then  compared  to  tables  of  quantiles  of  T  predicted  from 
a  postulated  T  distribution.  Table  B.l  shows  the  quantiles 
of  the  Kolmogorov  test  statistic.  The  quantity  n  is  the 
number  of  points  in  the  sample  and  p(=l-a)  is  the  "confi¬ 
dence  level". 

We  define  let  T^  be  the  test  statistic  associated 
with  one  individual  test.  The  first  quantile  value  higher 
than  Ti  is  associated  with  the  approximate  critical  con¬ 
fidence  and  significance  levels,  p  and  a.  For  example, 
if  the  value  of  Tj ,  based  on  a  sample  of  30  points,  were 
0.200,  then  0.8  <  $  <  0.9  and  0.2  >  £  >  0.1.  On  the  other 
hand,  if  Tj  for  the  same  case  were  0.288,  then  0.98  <  p  < 
0.99  and  0.02  >  ct  >  0.01.* 

Figure  4.1a-c  shows  the  graphical  comparison 
between  the  observed  and  the  calculated  X  ^  cumulative 
probabilities  for  YVETTE  08.  The  statistic  Tj  is  the 
greatest  separation  between  the  two  functions  along  a 
vertical  line  corresponding  to  one  of  the  observed  values  of 
S2.  in  all  the  cases  depicted  a  was  at  least  0.20.  Thus 
the  model  falls  within  the  80%  confidence  interval  of  the 
observed  distribution. 


*  The  quantity  p  *  1- a  defines  the  confidence  level  around 
the  empirical  distribution  at  which  the  model  distribution 
falls.  A  value  of  0.1  <  a<  0.2  implies  that  the  model  falls 
within  the  80  to  90%  confidence  intervals. 
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Only  four  stations  could  divided  into  three 
depth  regimes.  In  the  other  eleven  the  mixed  layer  was 
not  well  established.  In  some  cases,  the  layer  was  so  thin 
that  shear  values  calculated  over  8m  and/or  16m  produced 
only  4  or  5  values.  The  test  was  not  applied  to  samples 
with  fewer  than  6  elements.  Table  B.2  contains  the  a  values 
calculated  for  all  stations. 

A  breakdown  of  these  results  by  depth,  separation, 
and  a  range  (Table  B.3)  shows  that  over  75%  of  the  samples 
have  a  critical  level  of  0.20  or  higher.  Correspondingly, 
the  model  falls  within  the  80%  confidence  interval  for  75% 
of  the  samples. 

The  conservative  bias  of  the  goodness-of-f it 
test  could  have  been  avoided  by  using  the  x2  goodness- 
of-f  it  test  rather  than  the  Kolmogorov  test.  However,  that 
test  has  its  own  set  of  drawbacks.  For  one  thing,  the  x2 
test  assumes  a  sample  size  large  enough  that  the  test 
statistic  is  governed  by  a  x 2  distribution.  This  is  prob¬ 
ably  not  the  case  for  many  of  the  regimes  studied  here.  In 
contrast,  the  Kolmogorov  test  is  effective  for  small  sample 
size.  Secondly,  the  x2  test  depends  on  breaking  the 
sample  up  into  bins  of  arbitrary  width  and,  especially  with 
a  small  sample,  the  results  are  sensitive  to  bin  size. 
Lastly,  quantiles  of  the  x^  statistic  are  not  known  with  a 
great  deal  of  confidence  (Conover,  1980).  This  can  reduce 
the  usefulness  of  the  test.  For  these  reasons  we  chose  to 
accept  the  overestimated  values  predicted  by  the  Kolmogorov 
test  and  assume  that  the  errors  involved  are  lower  than 
those  inherent  in  the  use  of  the  X2  test  for  small  sample 
size. 
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4.2 


TESTING  MODEL  ASSUMPTIONS 


As  we  mentioned  in  Section  3,  the  reasoning  which 
led  to  the  hypothesis  of  a  y?  distribution  was  based  on 
physical  assumptions  about  the  nature  of  internal  wave- 
induced  shear.  We  can  apply  another  set  of  statistical 
tests  to  the  distributions  of  Azu  and  Azv  in  order  to 
test  how  well  these  satisfy  the  assumptions.  This  second 
approach  provides  an  independent  test  of  the  hypothesis. 

We  have  assumed  in  our  model  that  the  velocity 
difference  components  meet  the  following  assumptions: 

•  Azu  an(j  Azv  occur  randomly; 

•  they  are  distributed  normally; 


•  they  are  independent  quantities; 

•  they  come  from  identically  distributed  pop- 
uations . 

The  reasons  for  these  assumptions  were  discussed  in  Section 
3.  In  this  section  we  will  show  the  statistical  techniques 
we  used  to  determine  how  well  our  samples  satisfied  them. 

4.2.1  Cox  and  Stuart  Trend  Test 


We  did  not  attempt  to  assess  directly  the  random¬ 
ness  of  the  shear  component  samples.  Instead,  we  tested  to 
see  whether  or  not  a  trend  existed  using  the  Cox  and  Stuart 
Test.  A  sequence  of  numbers  will  have  a  trend  if  values 
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later  in  the  sequence  are  uniformly  higher  or  lower  than 
earlier  values.  If  no  such  tendency  exists,  then  the  sample 
is  probably  made  up  of  random,  independent  values. 

The  Cox  and  Stuart  test  consists  of  grouping  a 
sequence  of  variables,  of  length  n,  into  pairs  (x^,  x^+ij), 

where  k=n/2  if  n  is  even  and  k=(n  +  l)/2  if  n  is  odd.  The 

pairs  are  tested  to  see  which  value  is  larger.  If  we  asign 
a  "+”  to  pairs  in  which  x^  +  k  >  x^  and  a  to  pairs  in 

which  Xi+k  <  Xj ,  then  the  test  statistic  is  T,  the 

number  of  pairs  which  can  be  denoted  by  a  "  +  Our  null 

hypothesis  is  that  the  probability  of  observing  a  "+"  is  the 
same  as  the  probability  of  observing  a 

The  range  of  acceptance  is  bounded  according  to 
m-t  _<  T  _<  t ,  where  m  is  the  number  of  pairs  excluding  ties 

and  t  is  the  greatest  number  of  +'s  which  would  be  observed 

at  level  a  based  on  a’  probability  of  0.5.  (See  Conover 
(1980)  for  the  technique  for  calculating  t.)  If  T  falls 

outside  the  bounds,  then  the  hypothesis  is  rejected. 
The  results  of  this  test  are  shown  in  Tables  B.4  and  B.5. 
Tables  B.6  and  B.7  show  the  breakdown  by  regime,  separation, 
and  a.  A  total  of  91  Azu  profiles  and  98  Azv  profiles  (out 

of  114)  satisfied  the  test  at  S  >  0.10.  Somewhat  more  Azv 

profiles  than  Azu  profiles  satisified  at  a  levels  above  0.05, 
but  both  types  of  profiles  can  definitely  be  said  to  be  free 
of  trends. 


*  Ties,  where  x^+k  =  Xi,  are  disregarded  and  the  total 
number  of  pairs  is  reduced  by  the  number  of  ties. 
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Lilliefors  Normality  Test 


In  order  to  evaluate  whether  the  shear  com¬ 
ponent  samples  come  from  normal  distributions,  we  applied 
the  Lilliefors  test.  The  test  is  performed  not  on  an 
original  sequence  x^,  but  on  a  normalized  sequence. 

Zj.  —  ( *i  x  )  /  x  2  i  =  1 ,2, .  . .  ,n 

where  x^  is  the  sample  mean  and  X2  the  sample  variance. 
The  distribution  of  this  normalized  sequence  is  compared 
against  a  normal  distribution  with  zero  mean  and  unit 
variance.  If  we  define  F(z)  as  our  hypothesized  normal 
distribution  and  S(z)  as  the  observed  distribution  of  the 
adjusted  sequence,  then  the  test  statistic  is  the  same  as 
for  the  Kolmogorov  Test: 

suo , 

T  =  z  J  |F(z)  -  S(z) | . 

The  Kolomogorov  tables  cannot  be  used  to  evaluate  this 
test,  however;  instead,  special  quantiles  must  be  used  to 
make  a  decision. 

The  critical  levels  obtained  by  applying  this 
test  to  YVETTE  stations  are  shown  in  Tables  B.8,  and  B.9. 
The  distributions  of  a  values  are  listed  in  Tables  B.10  and 
B.ll.  The  test  results  showed  that  the  normality  assumption 
was  satisfied  by  nearly  all  profiles  (95  of  the  azu  and  96 
of  the  Azv)  at  levels  above  0.10. 


The  chi-squared  independence  test  assumes  that 
we  have  a  sample  of  length  N  whose  elements  can  be  sorted 
♦  by  two  criteria.  In  our  application  the  two  criteria  are 


whether  a  value  comes  from  A  zu  or  Azv  and  which  of  a 
number  of  bins  of  width  AS^  a  shear  component  value  fits. 
If  we  break  the  range  of  shear  into  C  bins,  then  each 
observation  can  only  fall  into  one  of  2C  categories. 

The  test  statistic  is  based  on  the  number  of 
entries  in  each  category.  We  can  envision  the  possible 
choices  as  a  matrix  with  2  rows  (one  for  each  shear  compo¬ 
nent)  and  C  columns  (one  for  each  shear  range  bin).  Each 
possible  category  can  be  denoted  by  indices  i  and  j,  where 
i=l  or  2,  j=l,2,....,  C.  The  number  of  observations  in 
each  bin  can  be  denoted  as  O^j.  For  each  row  we  define: 


L 


The  assumption  in  tnis  test  is  that  T  is  governed  by  a 
chi-square  distribution  with  (C-l)  degrees  of  freedom.* 
If  T  falls  below  the  a  level  quantile  then  we  can  accept  the 
hypothesis  that  an  element  occurs  in  row  1,  column  j 
independently  of  an  element  in  row  2,  column  j  and  that 
this  is  true  for  all  j.  For  our  case  this  is  equivalent  to 
saying  that  the  Au/Az  values  occur  independently  of  the 
Av/Az  values. 

The  critical  levels  obtained  by  applying  this 
test  to  YVETTE  profiles  are  listed  in  Table  B.12  and  the 
a  distribution  in  Table  B.13.  Nearly  all  of  the  profile 
pairs  (112)  satisfied  the  independence  test  at  a  levels 
above  0.10. 

This  test  does  not  seem  subject  to  the  bin- 
width  sensitivity  mentioned  when  discussing  the  chi-square 
goodness  of  fit  test.  We  presume  that  this  is  because  two 
observed  distributions  are  being  compared  with  one  another 
rather  than  one  observed  distribution  with  a  hypothesized 
distribution.  Nevertheless,  the  number  of  bins  was  never 
greater  than  one-half  the  number  of  samples,  and  in  most 
cases  was  less  than  one-third. 

4.2.4  The  Smirnov  Identical  Distribution  Test 


The  final  assumption  to  be  tested  is  that  the  two 
shear  components  are  governed  by  identical  distributions. 


*  Actually  with  (number  of  rows-1 ). (number  of  columns-1) 
degrees  of  freedom.  In  this  application,  with  only  two 
rows,  we  find  (C-l)  degrees  of  freedom. 
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The  test  chosen,  the  Smirnov  test,  is  of  the  Kolmogorov  - 
type  in  which  the  test  statistic  is  the  maximum  separation 
between  two  observed  distributions.  If  S^Cx)  and  S2(y)  are 
empirical  distributions  of  quantities  x  and  y,  and  Fj(x)  and 
F2(y)  are  their  (unknown)  hypothesized  distributions,  then 
we  can  test  the  hypothesis  F^(x)  =  F2(y)  by  generating  the 
statistic 


T  =S2P!  SX ( z>  -  S2  (z)l 

for  z's  in  the  range  of  all  possible  x’s  and  y's.  If  T 
is  less  than  the  appropriate  quantile,  then  the  hypothesis 

A 

can  be  accepted  at  level  a. 

Note  that  S^(x)  and  S2(y)  do  not  have  to  have 
the  same  number  of  samples  in  order  to  apply  this  test. 
However,  our  samples  are  of  equal  length. 

The  critical  levels  obtained  for  YVETTE  data 
are  listed  in  Table  B.14,  and  their  distribution  is  listed 
in  Table  B.15.  Again,  nearly  all  of  the  profile  pairs  (113) 
may  confidently  be  said  to  come  from  identical  distribu¬ 
tions. 


4.2.5  Summary  of  Assumption  Tests 

The  results  of  these  four  tests  show  that  we  are 
statistically  justified  in  accepting  the  hypotheses  based 
on  the  assumptions;  these  assumptions,  in  turn,  led  to  the 
postulated  x2  distribution  for  shear  squared.  That  is, 
we  can  assume  with  a  high  level  of  confidence  that  the  shear 
components,  Azu  and  Azv,  do  come  from  distributions  which 
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are  random,  normally  distributed,  mutually  independent,  and 
identically  distributed.  Table  4.1  shows  the  percentage  of 
each  set  of  samples  which  have  specific  minimum  critical 
levels.  For  example,  the  Cox-Stuart  Test  on  A  zu  was 
satisfied  by  69%  of  the  samples  at  a>0.20,  and  80%  at 
a>0.10.  These  results  lead  us  to  believe  that: 

•  the  agreement  between  S2  distribution  and 
X2  distributions  based  on  sample  variances  is 
significant,  and, 

•  that  because  the  assumption  tests  are  not 
biased,  we  could  extend  the  previous  statement 
to  indicate  that  S2  values  do  seem  to  come 
from  a  X2  distributed  population. 
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Table  4.1 

Percentages  of  Samples  which  satisfy 
assumption  tests  at  three  &  levels. 


Test 

a 

= 

0.20 

0.10 

0.05 

Cox-Stuart  Azu 

69% 

80% 

87% 

Cox-Stuart  Azv 
Lilliefors  Azu 
Lilliefors  Azv 
X2  Independence 


Smirnov 


Section  5 
MODEL  APPLICATION 


Thus  far  we  have  described  a  model  for  the 
distribution  of  shear  in  several  depth  regimes  of  the  upper 
ocean  and  established  a  statistical  basis  for  accepting  the 
model.  The  final  question  which  must  be  addressed  is  how  to 
apply  the  model . 

The  key  parameter  in  the  model  is  o^,  the  variance 
of  the  velocity  component  differences  based  on  observations. 
If  we  can  relate  this  parameter  to  large  scale,  more  easily 
observed  quantities,  then  we  can  develop  the  S2  probability 
distribution  that  we  seek.  One  parameter  which  can  be 
obtained  with  relative  ease  is  N2 ,  the  Brunt-V'ai s'ala 
frequency  at  depth  in  the  water  column.  Since  we  divided 
the  YVETTE  shear  profiles  into  regimes  distinguished  by 
changes  in  N2 ,  we  examined  the  possibility  of  obtaining 
values  of  through  some  correlation  between  N2  and  S2. 

Patterson  et  al.  (1981)  have  analyzed  the  same 
YVETTE  data  that  we  used  and  have  reached  some  conclusions 
about  the  relationship  between  N2  and  S2.  They  showed 
that  if  N2  and  S2  are  computed  over  2  m  differencing 
intervals,  then  the  ratio  of  N2  to  S2  (where  the  overbar 
indicates  depth  averaging  as  before)  varies  between  0.5 
and  3.3  in  the  upper  thermocline  and  deeper  layers  (Figure 
5.1) . 
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3.31 


U  -  Surface  Jlixed  Layer 
S  -  Near-Surface  Thennocline 
B  -  Belcw  Thennocline 
P  -  Subsurface  Pycnostad 
W  -  Whole  Profile 
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N1  M  11  10  12  25  8  24  9  5  21  23  18  EG  1J6 

YVETTE  Station  Number 


Figure  5.1 


~ 2  “5 

Ratio  of  N  to  S  (based  on  2  m  differencing 
intervals)  versus  station  number  for  the 
different  stratification  regimes  (from  Patterson 
et  ad. ,  1981) . 
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The  observation  that  this  ratio  is  reasonably 
well-behaved  motivated  us  to  estimate  as 

MS2  =  k  N2.  (5.1) 

where  k  is  an  order  unity  proportionality  constant.  We 
chose  a  value  of  k=l  as  a  first  approximation.  We  then 
estimated  the  sample  component  variance  from  (3.6b)  as 


We  tested  the  performance  of  the  x  2  model  based 
on  (5.2)  by  repeating  the  Kolmogorov  goodness-of-f it  test 
for  ten  of  the  YVETTE  stations.  The  values  of  the  test 
statistics,  T,  are  shown  in  Table  5.1,  along  with  the 
statistics,  Tq.201  which  is  the  maximum  T  value  which 
would  have  resulted  in  a  fit  at  level  of  significance 

A 

a  =  0.20.  If  we  had  been  testing  a  hypothesis  of  the  form 

Ho  =  "S^  occur  according  to  a  x2  -distribution 
with  variance  given  by  (5.2)," 

then  these  results  would  not  be  very  impressive.  Only  6  out 
of  20  profiles  produced  T  values  sufficient  to  result  in 
acceptance  at  a  comfortably  high  significance  level  of 

A 

a  2  0.20.  On  the  other  hand,  the  largest  calculated  T  in 
Table  5.1  is  0.284.  The  quantity  T  is  defined  as  the  maxi¬ 
mum  difference  between  the  observed  probability  distribution 
and  the  x2  distribution.  Thus  the  model-predicted  proba¬ 
bility  of  occurrence  (expressed  as  percentage  of  occurrence) 
of  S2  values  less  than  specified  levels  differed  from  the 


*  Table  5.1 

Test  Statistics  for  the  Kolmogorov  Test.  Test 
statistics,  T,  obtained  by  estimating  X 2 ( S f  og2 
based  on  (5.2).  The  calculations  were  only 
carried  out  for  ^z-2  m.  Tq.20  1S  the  statistic 
which  would  specify  rejection  at  level  0.20. 


STATION 

REGIME 

T 

YVETTE  05 

UT 

0.284 

DL 

0.114 

08 

UT 

0.065 

DL 

0.174 

09 

UT 

0.079 

DL 

0.114 

10 

UT 

0.062 

DL 

0.261 

11 

UT 

0.161 

DL 

0.183 

12 

UT 

0.178 

DL 

0.190 

18 

UT 

0.231 

DL 

0.182 

21 

UT 

0.183 

DL 

0.075 

23 

UT 

0.138 

DL 

0.145 

24 

UT 

0.260 

DL 

0.262 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


) .  20 

.151 

.081 

.138 

.096 

.102 

.102 

.126 

.084 

.130 

.084 

.139 

.079 

.111 

.091 

.098 

.103 

.145 

.078 

.187 
.296 
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observed  probability  by  no  more  than  30%  in  any  of  the 
profiles  tested.  And  in  75%  of  the  cases  the  difference 
was  less  than  20%.  These  figures  serve  to  indicate  the 
degree  of  confidence  that  the  user  can  place  in  predictions 
from  the  statistical  model  used  in  this  manner. 


Section  6 
SUMMARY 


We  have  described  a  statistical  model  for  the 
distribution  of  S2  values  in  different  depth  regimes  in 
the  upper  ocean.  The  model  indicates  that  a  distribution 
with  variance  computed  from  the  average  S2  in  the  sample 
is  a  good  approximation  to  the  observed  distribution.  The 
model  is  physically  reasonable  and  has  been  verified  by 
statistical  analysis  of  shear  profiles  measured  at  15 
different  YVETTE  stations. 

The  key  to  translating  this  model  into  a  useful 
tool  for  future  study  lies  in  developing  a  method  of 
estimating  the  mean  value  of  S2 ,  PS2,  in  a  depth  regime 
from  easily  measured  oceanic  parameters.  A  simple  test, 
based  on  results  of  an  independent  analysis  of  YVETTE  data, 
consisted  of  estimating  ^S2  from  the  mean  Brunt- Vaisala 
frequency  in  the  regime.  This  test  showed  that  such  a 
simplified  approach  to  using  the  statistical  model  could 
yield  a  cumulative  shear  distribution  which  agreed  with  the 
observed  distribution  to  30%  or  better.  This  particular  use 
of  the  model  is  restricted  at  present  to  shear  values 
obtained  from  differences  over  two  meters. 

The  reasonable  success  of  the  X2  model  when 
combined  with  a  fairly  primitive  scheme  for  estimating 
the  parameter  needed  by  the  model  leads  us  to  suggest 
further  study.  This  would  include  examination  of  more  S2 
data  sets,  both  vertical  profiles  and  time  series  from 
vertically  spaced  current  meters.  We  would  also  need 
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rtirno^i 


to  examine  more  extensively  the  nature  of  the  N2-S2  rela¬ 
tionship  over  various  vertical  scales.  In  particular, 
it  appears  (Patterson  et  al . ,  1981)  that  the  relationship 

between  n2  and  S2  may  be  less  well  defined  at  larger 
separations.  This  would  of  course,  reduce  the  effective¬ 
ness  of  the  model  over  these  scales. 


APPENDIX  A 
HYPOTHESIS  TESTING 


It  is  possible  to  employ  statistical  tests  to 
quantify  the  degree  of  confidence  which  can  be  assigned  to  a 
hypothesis  such  as  ”S2  comes  from  a  distribution".  To 
do  this  we  establish  a  null  hypothesis,  H0 ,  regarding  some 
characteristic  of  the  data  set  being  examined: 

H0  =  "the  hypothesis  being  tested  is  true." 

We  then  use  the  data  to  produce  a  test  statistic,  T,  whose 
value  will  allow  us  to  either  reject  or  not  reject  H0. 
Notice  the  choice  of  words  in  the  last  sentence.  If  we 
reject  H0 ,  we  are  in  effect  stating  "the  hypothesis  is 
false."  But  if  we  do  not  reject  H0  we  are  not  saying 
"the  hypothesis  is  true."  Instead,  we  are  saying  "the 
hypothesis  cannot  be  proven  false."  •  We  must  allow  for  the 
possibility  that  a  larger  sample  might  lead  us  to  reject 
H0 .  Tests  of  this  type  are  rejection  tests  in  that  we 
can  only  reject  or  not  reject,  but  never  accept  H0  without 
reservation.  It  is  conceivable,  for  example,  that,  based 
on  the  analysis  of  a  single  data  set,  two  mutually 
exclusive  hypotheses  could  both  fail  to  be  rejected.  A 
detailed  description  of  hypothesis  testing  techniques 
is  presented  by  Conover  (1980).  We  present  an  outline 
here. 

Hypothesis  testing  is  subject  to  two  types  of 

error : 

•  rejecting  H0  when  H0  is  actually  true  (Type  I 
error) ,  and 
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•  accepting  H0  when  H0  is  actually  false  (Type  II 
error) . 

We  can  quantify  approximate  probabilities  of 
making  these  errors  based  on  a  given  sample.  We  define  a  , 
the  level  of  significance,  as  the  maximum  probability  of 
making  a  Type  I  error,  and  6  as  the  probability  of  making  a 
Type  II  error. 

We  want  both  a  and  S  to  be  as  small  as  possible. 
Unfortunately,  for  a  fixed  sample  size,  decreases  in  a 
are  usually  accompanied  by  undesirable  increases  in  B  . 
That  is,  if  we  set  up  a  test  which  tends  to  minimize  the 
probability  of  rejecting  a  true  null  hypothesis,  we  increase 
the  chance  of  accepting  a  false  null  hypothesis.  The  only 
way  to  reduce  both  a  and  B  is  to  increase  the  sample  size 
(until  eventually  the  sample  encompasses  all  members  of  the 
population  being  studied  and  all  tests  are  exact).  Only  if 
there  is  a  simple  alternative  hypothesis  to  H0  can  £  be 
estimated.  This  is  not  the  case  in  general. 

The  rationale  behind  the  level  of  significance 
can  be  illustrated  by  a  simple  example.  Suppose  we  have 
a  finite  sample  of  n  random  independent  observations  of 
some  quantity  x  from  which  we  compute  a  given  estimator  T. 
We  have  reason  to  believe  that  the  true  value  T  (given 
infinitely  many  observations)  of  the  parameter  being  esti¬ 
mated  is  T0 .  Thus  our  hypothesis  H0  is  T=T0.  Even  if 
T=T0,  the  sample  value  T  will  not  equal  T0 .  We  assume 
that  a  sample  value  T  has  a  probability  density  distribution 
B(T) ;  then  we  can  use  this  distribution  to  estimate  the 
pr  -'ability  of  observing  any  value  of  T  given  that  H0  is 
true . 
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Figure  A-l  shows  a  schematic  probability  density 
which  assumes  T=T0 .  The  test  statistic  from  a  single 
test,  Ti ,  falls  some  distance  from  T0 .  If  is  so  far 
from  T0  that  its  probability  of  occurrence  is  very  small, 
we  decide  that  H0  should  be  rejected.  In  order  to  establish 
criteria  for  our  decision  we  decide  on  upper  and  lower 


limits  —a/ 2  and  -1-0/2' 


which  define  a  region  of  acceptance. 


These  T^q  are  chosen  such  that  the  probabilities  of  values 
from  the  distribution  falling  beyong  these  limits  are 


Prob  (T  <  T1_^a)  =  $a 


(A. la) 


and 

Prob  (T  >  T4a)  =  Ja. 


(A. lb) 


Thus  the  probability  of  T  falling  outside  the  range  defined 
by  these  limits  is  a.  We  set  a  small;  then  if  we  do  observe 
such  a  T  we  can  be  confident  that  we  have  only  a  low  prob¬ 
ability  of  being  wrong  by  rejecting  H0.  a  is  called 
the  "level  of  significance"  of  the  test  and  the  range  of  T 
for  which  the  hypothesis  will  be  rejected  is  called  the 
"region  of  rejection."  If  H0  is  true,  the  chance  that  Tj_ 
will  fall "in  the  range  of  rejection  is  just  a;  thus  a  is  the 
probability  of  rejecting  a  true  null  hypothesis  (Type  I 
error) . 


We  can  also  define  a  second  quantity,  a,  which 
we  refer  to  as  the  critical  level  of  significance.  This 
is  defined  as  the  smallest  significance  level  at  which 
the  hypothesis  would  be  rejected  for  a  given  observation. 
This  is  defined  by 


Prob  (T  >  Tx )  =  a/2 


9 


(A. 2) 
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Figure  A.l 


Schematic  probability  density  of  T.  The 
double  cross-hatching  indicates  the  area 
of  rejection  defined  by  a  .  The  single- 
and  double-cross-hatching  together  indi¬ 
cate  the  critical  area  defined  by  <£. 
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if  >  T0,  and  a  similar  definition  if  Tj  <  T0.  This  level 

is  also  illustrated  in  the  figure.  If,  as  illustrated, 
>\  « 
a  >  a,  we  have  the  situation  in  which  the  level  of  signifi¬ 
cance  criterion  is  more  than  satisfied.  Thus  a  can  be 
looked  at  as  a  measure  of  the  robustness  of  the  test. 
Asa  increases,  the  confidence  with  which  we  can  decide  not 
to  reject  H0  increases. 

It  is  important  to  notice  that  a  is  arbitrarily 
defined  by  the  tester,  while  a  is  defined  by  the  sample 
being  tested.  We  can  require  that  the  test  be  satisfied  at 
some  small  level  a.  If  we  do  this,  we  are  taking  the  point 
of  view  that  the  sample  must  have  every  reasonable  opportun¬ 
ity  to  pass  the  test  and  that  we  want  to  minimize  our 
chances  of  rejecting  a  true  H0.  Or  we  can  let  the  sample 
indicate  how  well  or  poorly  it  satisifes  the  test  by  means 
of  a  .  This  way,  our  point  of  view  is  that  we  are  less 
concerned  with  the  possibility  of  rejecting  a  true  hypoth¬ 
esis  than  we  are  with  assessing  the  certainty  with  which  we 

/\ 

make  decisions.  We  choose  to  use  a  in  our  analysis. 


t 
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APPENDIX  B 

STATISTICAL  RESULTS  BASED  ON 
TESTS  OF  YVETTE  DATA 


Table  B.l 


Quantiles  of  the  Kolmogorov  Test  Statistic*. 


Two-SM'd  Ten 


•90 

M 

.91 

•99 

J 

•90 

-95 

•99 

S9 

m  rnm  1 

.900 

.950 

.975 

.990 

.995  ii  - 

21 

.226 

.259 

.287 

.321 

.344 

2 

.644 

.776 

.842 

.900 

.929 

22 

.221 

.253 

.281 

.314 

.337 

3 

.565 

.<36 

.708 

.785 

.829 

23 

.216 

J47 

.275 

.307 

.330 

4 

.493 

.565 

.624 

.689 

.734 

24 

.212 

.242 

J69 

.301 

.323 

3 

.447 

.509 

.563 

.627  * 

.669 

25 

JOS 

J38 

.264 

J95 

.317 

6 

.410 

.466 

.519 

.577 

.617 

26 

.204 

.233 

J59 

.290 

.311 

7 

.3S1 

.436 

.483 

.538 

.576 

27 

.200 

.229 

.254 

.284 

.305 

S 

.351 

.410 

.454 

.507 

342 

28 

.197 

.225 

.250 

.279 

.300 

9 

.339 

.387 

.430 

.480 

.513 

29 

.193 

.221 

.246 

.275 

.295 

10 

.323 

.369 

.409 

.457 

.489 

30 

.190 

.218 

.242 

.270 

.290 

11 

.301 

.352 

.391 

.437 

.468 

31 

.187 

.214 

.238 

.266 

.285 

12 

.296 

.338 

.375 

.419 

.449 

32 

.184 

.211 

.234 

.262 

.281 

13 

.255 

.325,. 

.404 

.432 

33 

.182 

.201. 

.231 

.258 

J77 

14 

.275 

.3147 

.349 

.390 

.418  ' 

34 

.179 

.205 

.227 

.254 

.273 

15 

.266 

.304 

.338 

.377 

.404 

35 

.177 

.202 

.224 

.251 

J69 

16 

.251  * 

.295  • 

:.’327 

.366 

.392 

36 

.174 

.199 

.22! 

.247 

J6S 

17 

.250 

.286 

.318 

.355 

.381 

37 

.172 

.196 

J18 

.244 

.262 

It 

.244 

.279 

.309 

.346 

.371 

38 

.170 

.194 

.215 

.241 

J58 

19 

.237 

.271 

.301 

.337 

.361 

39 

.168 

.191 

.213 

.238 

.255 

20 

.232 

.265 

.294 

.329 

.352 

40 

.165 

.189 

.210 

J35 

.252 

Approximation 
for  n  >  40 

1.07 

vl 

1.22 

vS 

1.36 

"vl 

1.52 

v  J* 

1.63 

v I 

The  entries  in  this  table  are  selected  quantiles 
of  the  Kolmogorov  test  statistics  T.  Reject  H0  at  the 
level  a  if  T  exceeds  the  1-a  quantile  given  in  this  table. 
These  quantiles  are  exact  for  n  £  40.  The  other  quantiles 
are  approximations  that  are  equal  to  the  exact  quantiles  in 
most  cases.  A  better  approximation  for  n  >  40  results  if 
(n+  V  n/10 ) ^  is  used  instead  of  VtT  in  the  denominator. 
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Table  B.2 

Kolmogorov  Good ness-of -Fit  Test  Results 


Critical  levels*  for  the  three  depth  regimes  for  all  15  YVETTE  stations. 


Depth  Regime  =  ML  UT  DL 


Separation  = 

’ 

2m 

4 

8 

16 

2 

4 

8 

16 

2 

4 

8 

16 

Station 

05  . 

<0.01 

0.20 

— 

— 

<0.01 

0.02 

0.20 

0.20 

0.20 

0.10 

0.20 

08 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

09 

0.20 

0.20 

0.01 

— 

0.20 

0.10 

0.10 

0.20 

0.20 

0.20 

0.20 

0.20 

10 

— 

— 

— 

— 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

0.05 

0.10 

11 

— 

— 

— 

-  • 

0.02 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

12 

— 

— 

— 

— 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

18 

— 

— 

— 

— 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

21 

— 

— 

— 

— 

0.02 

0.05 

0.20 

<0.01 

0.20 

0.20 

0.20 

23 

— 

— 

— 

— 

0.20 

0.20 

0.20 

<0.01 

<0.01 

<0.01 

0.02 

24 

— 

— 

— 

— 

0.10 

0.20 

— 

— 

0.20 

0.05 

— 

— 

25 

— 

— 

— 

— 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

N0R1 

0.20 

0.20 

— 

— 

0.05 

0.20 

0.02 

— 

N0R4 

— 

— 

— 

— 

0.20 

— 

— 

— 

0.20 

0.20 

0.20 

— 

NOR6 

— 

— 

— 

— 

<0.01 

<0.01 

<0.05 

0.20 

— 

— 

— 

— 

EP0CS6 

0.10 

0.20 

— 

— 

0.20 

0.20 

0.20 

0.20 

0.02 

0.20 

0.20 

0.20 

4 

4 

2 

1 

15 

13 

12 

12 

14 

14 

13 

11 

^  — 

N— 

■ 

> 

w 

11  52  52 


*  &  values  listed  are  not  exact.  Table  A.l  shows  five 
levels:  0.20,  0.10,  0.05,  0.02,  0.01.  The  &  value  listed 

in  this  table  indicates  the  approximate  range  of  a.  For 
example,  a  value  of  0.10  implies  0.10  <&<  0.20. 
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Table  B.3 

Kolmogorov  Goodness-of-Fit  Test  Results 


Distribution  of  critical  levels  by  a  range  for  depth 
regimes  and  separations  from  all  YVETTE  stations. 


a  >0.2 

o.2>a>o.i 

o.i>a_>o.o5 

0.05>a>0.02 

o.o2>a>o.oi 

a <0.01 

TOTAI 

Az  = 

2m 

2 

i 

0 

0 

0 

i 

4 

ML 

= 

4m 

4 

0 

0 

0 

0 

0 

4 

= 

8m 

1 

0 

0 

0 

1 

0 

2 

= 

16m 

1 

0 

0 

0 

0 

0 

1 

A  z  = 

2m 

10 

1 

0 

1 

1 

2 

15 

ITT 

= 

4m 

10 

1 

0 

1 

1 

0 

13 

= 

8m 

9 

1 

2 

0 

0 

0 

12 

s 

16m 

12 

0 

0 

0 

0 

0 

12 

Table  B.4 

Cox-Stuart  Trend  Test  Results 


Critical  levels  for  the  Azu  distri¬ 
butions  for  all  YVETTE  stations. 


DEPTH  REGIME  =  ML 


UT 


DL 


SEPARATION  = 

=  2 

4 

8 

16 

2 

4 

8 

16 

2 

4 

8 

16 

YVETTE  05 

0,20 

0.10 

— 

— 

0.20 

0.20 

0.10 

— 

0.20 

0.20 

0.20 

0.20 

08 

0.20 

0.01 

0.20 

0.10 

<0.01 

0.05 

0.01 

0.05 

0.20 

0.20 

0.20 

0.10 

09 

0.20 

0.05 

0.20 

— 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

10 

— 

— 

— 

— 

0.01 

<0.01 

0.10 

0.05 

<0.01 

<0.01 

0.02 

0.02 

11 

— 

— 

— 

— 

<0.01 

0.05 

0.05 

0.02 

0.02 

0.20 

0.10 

0.20 

12 

— 

— 

— 

— 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

18 

— 

— 

— 

— 

0.20 

0.20 

0.20 

0.20 

0.02 

0.20 

0.20 

0.20 

21 

— 

— 

— 

— 

<0.01 

0.20 

0.10 

0.20 

0.10 

0.20 

0.20 

0.20 

23 

— 

— 

— 

— 

0.20 

0.05 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

24 

— 

— 

— 

— 

0.10 

0.20 

— 

— 

0.20 

0.20 

— 

— 

25 

— 

— 

— 

— 

0.20 

0.20 

0.20 

0.10 

0.20 

0.20 

0.20 

0.20 

NOR 

1 

— 

— 

— 

— 

0.20 

— 

— 

0.20 

0.20 

0.20 

— 

4 

— 

— 

— 

— 

0.10 

— 

— 

— 

0.20 

0.20 

— 

— 

6 

— 

— 

— 

— 

0.20 

0.05 

0.20 

0.20 

— 

— 

— 

— 

6 

0.20 

0.10 

— 

— 

0.20 

0.20 

0.20 

0.20 

0.01 

0.20 

0.20 

0.20 

B-4 


Table  B.5 

Cox-Stuart  Trend  Test  Results 


Critical  levels  for  the  Azv  distri¬ 
butions  for  all  YVETTE 


DEPTH  REGIME  =  ML 


(JT 


DL 


SEPARATION  = 

2 

4 

8 

16 

2 

4 

8 

16 

2 

4 

8 

16 

YVETTE 

05 

0.20 

0.10 

_ 

_ 

0.20 

0.20 

0.05 

__ 

0.20 

0.20 

0.20 

0.20 

08 

0.05 

0.20 

0.10 

0.20 

0.20 

0.20 

0.20 

0.05 

0.20 

0.20 

0.20 

0.10 

09 

0.20 

0.20 

0.20 

— 

0.20 

0.20 

0.10 

0.20 

0.20 

0.20 

0.20 

0.20 

10 

— 

— 

— 

| 

0.05 

0.20 

0.10 

0.20 

0.20 

0.20 

0.20 

0.20 

11 

— 

— 

— 

— 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

12 

— 

— 

— 

— 

0.10 

0.20 

0.20 

0.02 

<0.01 

0.10 

0.10 

0.05 

18 

— 

— 

— 

— 

0.20 

0.20 

0.20 

0.10 

0.20 

0.20 

0.20 

0.20 

21 

— 

— 

— 

— 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

23 

— 

— 

— 

— 

0.20 

0.20 

0.20 

0.02 

0.20 

0.20 

0.20 

0.20 

24 

— 

— 

— 

— 

0.20 

0.20 

— 

— 

0.20 

0.02 

— 

— 

25 

— 

— 

— 

— 

0.20 

0.05 

0.20 

0.10 

0.20 

0.20 

0.20 

0.20 

NOR 

1 

— 

— 

— 

— 

0.20 

0.02 

— 

— 

0.02 

0.20 

0.20 

— 

4 

— 

— 

— 

— 

0.10 

— 

— 

— 

0.20 

0.20 

— 

— 

6 

— 

— 

— 

— 

0.20 

0.20 

0.20 

0.20 

— 

— 

— 

— 

EPOCS 

6 

0.20 

0.10 

— 

— 

<0.01 

0.05 

0.05 

0.02 

0.20 

0.20 

0.20 

0.20 

* 
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Table  B.6 

Cox-Stuart  Trend  Test  Results 


Distribution  of  critical  levels  for  Azu  by  fit  range  for  depth 
regimes  and  separations  for  all  YVETTE  stations. 


a>o.2o 

o.2o>a>o.io 

o.io>a>p.o5 

0.05>S>0.02 

0.02>S>0.01 

fit  <0.01 

TOTAL 

Az=2 

4 

0 

0 

0 

0 

0 

4 

4 

0 

2 

i 

0 

1 

0 

4 

ML  8 

2 

0 

0 

0 

0 

0 

2 

16 

0 

1 

0 

0 

0 

0 

1 

Az=2 

9 

2 

0 

0 

1 

3 

15 

UT  4 

9 

0 

4 

0 

0 

1 

14 

8 

7 

3 

1 

0 

1 

0 

12 

16 

7 

1 

2 

1 

0 

0 

11 

A  z=2 

9 

1 

0 

2 

1 

1 

14 

DL  4 

13 

0 

0 

0 

0 

1 

14 

8 

10 

1 

0 

1 

0 

0 

12 

16 

9 

1 

0 

1 

0 

0 

11 

TOTAL 

79 

12 

8 

5 

4 

6 

114 
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Table  b.7 

Cox-Stuart  Trend  Test  Results 

l 


Distribution 

of  critical 

levels 

for  Azu  by 

a  range  for 

depth 

regimes 

and  separations  for 

all  YVETTE 

stations . 

1 

a >0.20  0.20>S>0.10  0. 

10>a>0.05 

0.05>S>0.02 

0.02)6  >0.01 

S<0.01 

TOTAL 

1 

A  z=2 

3 

0 

1 

0 

0 

0 

4 

4 

2 

2 

0 

0 

0 

0 

4 

ML  8 

1 

1 

0 

0 

0 

0 

2 

16 

1 

0 

0 

0 

0 

0 

1 

> 

A  z=2 

11 

2 

1 

0 

_  0 

1 

15 

UT  4 

11 

0 

2 

1 

0 

0 

14 

8 

8 

2 

2 

0 

0 

0 

12 

» 

16 

5 

2 

1 

3 

0 

0 

11 

Az=2 

12 

0 

0 

1 

0 

1 

14 

DL  4 

12 

1 

0 

1 

0 

0 

14 

1 

8 

11 

1 

0 

0 

0 

0 

12 

16 

9 

1 

1 

0 

0 

0 

11 

TOTAL 

86 

12 

8 

6 

0 

2 

114 

I 
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Table  B.8 


Lilliefors  Normality  Test  Results 
Critical  levels  for  Azu  profiles  from  all  YVETTE  stations. 


DEPTH  REGIME  =  ML 


UT 


DL 


SEPARATION  = 

2 

4 

8 

16 

2 

4 

8 

16 

2 

4 

8 

16 

YVETTE 

:  05 

0.01 

<0.01 

~ 

j 

0.01 

0.01 

0,20 

—  i 

0.20 

0.20 

0.20 

0.20 

08 

0.05 

0.01 

0.20 

0.10 

0.05 

0.20 

0.01 

0.05 

0.20 

0.20 

0.20 

0.20 

09 

0.20 

0.10 

0.10 

1 

0.05 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

10 

— 

— 

— 

— 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

11 

— 

— 

— 

— 

0.10 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

12 

'  — 

— 

— 

— 

0.05 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

0.10 

18 

— 

— 

— 

— 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

21 

— 

— 

— 

— 

0.05 

0.05 

0.20 

0.20 

0.20 

0.20 

0.05 

0.20 

23 

— 

— 

— 

— 

0.20 

0.20 

0.20 

0.20 

<0.01 

0.10 

0.05 

0.20 

24 

— 

— 

— 

— 

0.10 

0.20 

— 

— 

0.20 

0.20 

— 

— 

25 

— 

— 

— 

— 

0.05 

0.20 

0.05 

0.20 

0.20 

0.20 

0.20 

0.20 

NOR 

1 

— 

— 

— 

— 

0.20 

0.20 

— 

— 

0.20 

0.20 

0.20 

— 

4 

— 

— 

— 

— 

0.20 

— 

— 

— 

0.20 

0.20 

— 

— 

6 

— 

— 

— 

— 

0.20 

0.20 

0.10 

0.20 

— 

— 

— 

— 

6 

0.20 

0.05 

— 

— 

0.20 

0.20 

0.20 

0.20 

0.10 

0.20 

0.20 

0.20 

Table  B.9 

Lilliefors  Normality  Test  Results 
Critical  levels  for  A zu  profiles  from  all  YVETTE  stations. 


DEPTH  REGIMES  =  ML  UT  DL 


SEPARATION  =248  16  248  16  24  8  16 


YVETTE 

I  05 

0.01 

0.05 

_ 

0.02 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

08 

0.20 

0.20 

0.20 

0.10 

0.20 

0.20 

0.20 

0.01 

0.20 

0.20 

0.20 

0.01 

09 

0.20 

0.20 

0.10 

— 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

10 

— 

— 

— 

— 

0.20 

0.20 

0.20 

0.10 

0.20 

0.20 

<0.01 

0.20 

11 

— 

— 

— 

— 

0.01 

0.05 

0.20 

0.05 

0.20 

0.20 

0.20 

0.20 

12 

— 

— 

— 

— 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

18 

— 

— 

— 

— 

0.20 

0.20 

0.20 

0.05 

0.20 

0.20 

0.20 

0.20 

21 

— 

— 

— 

— 

0.01 

0.20 

0.20 

0.20 

0.05 

0.20 

0.20 

0.20 

23 

— 

— 

— 

— 

0.20 

0.20 

0.20 

0.20 

<0.01 

0.10 

0.20 

0.20 

24 

— 

— 

— 

— 

0.20 

0.20 

— 

— 

0.20 

0.20 

— 

— 

25 

— 

— 

— 

— 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

NOR 

1 

— 

— 

— 

— 

0.10 

0.20 

— 

— 

0.20 

0.20 

0.20 

— 

4 

— 

— 

— 

— 

0.20 

— 

— 

— 

0.20 

0.20 

— 

— 

6 

— 

— 

— 

— 

<0.01 

0.10 

0.05 

0.20 

— 

— 

— 

— 

EPOCS 

6 

0.20 

0.20 

— 

— 

<0.05 

0.20 

0.10 

0.20 

0.20 

0.20 

0.20 

0.20 

Table  B.10 

Lilliefors  Normality  Test  Results 

Distribution  of  critical  levels  for  Azu  by  a  range  for 
depth  regimes  and  separations  for  all  YVETTE  stations. 


a  >0.20 

o.2o>a_>o.io 

o.io>a>o.o5 

0.05>a>0.02 

0.02>s>0.01 

a<o.oi 

TOTAL 

Az=2 

2 

0 

i 

0 

1 

0 

4 

4 

0 

i 

i 

0 

1 

1 

4 

ML  8 

1 

i 

0 

0 

0 

0 

2 

16 

0 

i 

0 

0 

0 

0 

1 

Az=2 

7 

2 

5 

0 

1 

0 

15 

UT  4 

12 

0 

1 

0 

1 

0 

14 

8 

9 

1 

1 

0 

1 

0 

12 

16 

10 

0 

1 

0 

0 

0 

11 

Az=2 

12 

1 

0 

0 

0 

1 

14 

DL  4 

13 

1 

0 

0 

0 

0 

14 

8 

10 

0 

2 

0 

0 

0 

12 

16 

10 

1 

0 

0 

0 

0 

11 

TOTAL 

86 

9 

12 

0 

5 

2 

114 
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Table  B.ll 


Lilliefors  Normality  Test  Results 

Distribution  of  critical  levels  for  A  zv  by  Q.  range  for 
depth  regimes  and  separations  for  all  YVETTE  stations. 


a>0 . 20 

o.2o>a>o.io 

o.io>a>o.o5 

0.05>a>0.02 

o.o2>a>p.oi 

a<o.oi 

TOTAL 

Az=2 

3 

0 

0 

0 

1 

0 

4 

4 

3 

0 

i 

0 

0 

0 

4 

ML  8 

1 

0 

i 

0 

0 

0 

2 

16 

0 

0 

i 

0 

0 

0 

1 

Az=2 

9 

i 

i 

1 

2 

1 

15 

UT  4 

12 

i 

i 

0 

0 

0 

14 

8 

10 

i 

i 

0 

0 

0 

12 

•  16 

7 

i 

2 

0 

1 

0 

11 

Az=2 

12 

0 

1 

0 

0 

1 

14 

DL  4 

13 

1 

0 

0 

0 

0 

14 

8 

11 

0 

0 

0 

0 

1 

12 

16 

10 

0 

0 

0 

1 

0 

11 

TOTAL 

91 

5 

9 

1 

5 

3 

114 

TOTAL 


Table  B.12 


X ^-Independence  Test  Results 
Critical  levels  for  profiles  from  all  YVETTE  stations. 


DEPTH  REGIME  =  ML 


UT 


DL 


SEPARATION  = 

2 

4 

8 

16 

2 

4 

8 

16 

2 

4 

8 

16 

YVETTE 

05 

0.10 

0.25 

- - 

0.25 

0.25 

0.25 

_ 

0.25 

0.25 

0.10 

0.25 

08 

0.25 

0.25 

0.25 

0.25 

0.25 

0.25 

0.25 

0.25 

0.25 

0.25 

0.25 

0.25 

09 

0.25 

0.10 

0.25 

— 

0.25 

0.25 

0.25 

0.25 

0.25 

0.25 

0.25 

0.25 

10 

— 

— 

— 

— 

0.25 

0.25 

0.25 

0.25 

0.25 

0.25  ' 

0.25 

0.25 

11 

— 

— 

— 

— 

0.25 

0.25 

0.10 

0.25 

0.25 

0.25 

0.25 

0.25 

12 

— 

— 

— 

— 

0.25 

0.25 

0.25 

0.25 

0.25 

0.25 

0.25 

0.25 

18 

— 

— 

— 

— 

0.25 

0.25 

0.25 

0.25 

0.25 

0.25 

0.25 

0.25 

21 

— 

— 

— 

— 

0.25 

0.25 

0.25 

0.25 

0.25 

0.25 

0.25 

0.25 

23 

— 

— 

— 

— 

0.25 

0.25 

0.25 

0.25 

0.25 

0.25 

0.25 

0.25 

24 

— 

— 

— 

— 

0.25 

0.25 

— 

— 

0.25 

0.25 

— 

25 

— 

— 

— 

— 

0.10 

0.25 

0.25 

0.25 

0.05 

0.10 

0.25 

0.25 

NOR 

1 

— 

— 

— 

— 

0.25 

0.25 

— 

— 

0.25 

0.25 

0.25 

— 

4 

— 

— 

— 

— 

0.25 

— 

— 

— 

0.02 

0.25 

— 

— 

6 

— 

— 

— 

— 

0.25 

0.25 

0.25 

0.25 

— 

— 

— 

— 

6 

0.25 

0.25 

— 

— 

0.25 

0.25 

0.25 

0.25 

0.25 

0.25 

0.25 

0.25 
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Table  B.13 

X^-Independence  Test  Results 

II 

Distribution  of  critical  levels  by  a  range  for  depth 
regimes  and  separations  for  all  YVETTE  stations. 


I 


3  >0.20 

0.20>q>0.10 

o.io>a>o.o5 

0.05>a>0.02 

0.02)a>p.01 

a  <0.01 

TOTAL 

Az=2 

3 

1 

0 

0 

0 

0 

4 

» 

4 

3 

1 

0 

0 

0 

0 

4 

ML  8 

2 

0 

0 

0 

0 

0 

2 

16 

1 

0 

0 

0 

0 

0 

1 

A  z=2 

14 

1 

0 

0 

0 

0 

15 

UT  4 

14 

0 

0 

0 

0 

0 

14 

8 

11 

1 

0 

0 

0 

0 

12 

16 

11 

0 

0 

0 

0 

0 

11 

• 

Az=2 

12 

0 

i 

1 

0 

0 

14 

DL  4 

13 

1 

0 

1 

0 

0 

14 

8 

11 

1 

0 

0 

0 

0 

12 

# 

16 

11 

0 

0 

0 

0 

0 

11 

TOTAL 

106 

6 

1 

1 

0 

0 
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Table  B .  14 

Smirnov  Identical  Distributions  Test  Results 


Critical  levels  for  profiles  from  all  YVETTE  stations 


0.20  0.20  0.20  — 

0.20  0.20  0.10  0.10 

0.20  0.20  0.20  0.20 

0.20  0.20  0.20  0.20 

0.10  0.20  0.20  0.20 

0.20  0.20  0.20  0.20 

0.20  0.20  0.20  0.20 

0.20  0.20  0.20  0.20 

0.20  0.20  0.20  0.20 

0.20  0.20  — 

0.20  0.20  0.20  0.20 

0.20  0.20  — 

0.20  —  — 

0.20  0.20  0.10  0.20 

0.20  0.20  0.20  0.20 


0.20  0.20  0.20  0.20 

0.20  0.20  0.20  0.20 

0.20  0.20  0.20  0.20 

0.20  0.20  0.20  0.20 

0.20  0.20  0.20  0.20 

0.20  0.20  0.20  0.20 

0.20  0.20  0.20  0.20 

0.20  0.20  0.20  0.20 

0.20  0.20  0.20  0.20 

0.20  0.20  — 

0.20  0.20  0.20  0.20 

0.20  0.20  0.20  — 

0.20  0.20  — 


) 


> 
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Table  B.15 

Smirnov  Identical  Distributions  Test  Results 

Distribution  of  &  for  depth  regimes  and  separations 
from  all  YVETTE  stations. 


> 

a  >0.20  0.20>S>0.10  0.10>i>0.05  0.05>6t>0.02  0.02>a>0.01  a  <0.01  TOTAL 


Ajs=2  3  0 

4  4  0 

ML  8  2  0 

16  1  0 


1 

0 

0 

0 


0 

0 

0 

0 


0  0  4 

0  0  4 

0  0  2 

0  0  1 


Az=2 

14 

1 

UT  4 

14 

0 

8 

10 

2 

16 

10 

1 

Az=2 

14 

0 

DL  4 

14 

0 

8 

12 

0 

16 

11 

0 

TOTAL 

109 

4 

.0 

0 

0 

0 

0 

0 

0 

0 

1 


0 

0 

0 

0 

0 

0 

0 

0 

0 


0  0  15 

0  0  14 

0  0  12 

0  0  11 

0  0  14 

0  0  14 

0  0  12 

0  0  11 

0  0  114 
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